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Digital  images  of  overcast  skies  as  seen  from  the  earth’s  surface  open  new  windows  onto  the  angular 
details  of  overcast  colors  and  visible-wavelength  spectra.  After  calibration  with  a  spectroradiometer,  a 
commercial  CCD  camera  equipped  with  a  fisheye  lens  can  produce  colorimetrically  accurate  all-sky  maps 
of  overcast  spectra.  Histograms  and  azimuthally  averaged  curves  of  the  resulting  chromaticities  show 
consistent,  but  unexpected,  patterns  in  time-averaged  overcast  colors.  Although  widely  used  models  such 
as  LOWTRAN7  and  MODTRAN4  cannot  explain  these  characteristic  patterns,  a  simple  semiempirical 
model  based  on  the  radiative  transfer  equation  does,  and  it  provides  insights  into  the  visible  consequences 
of  absorption  and  scattering  both  within  and  beneath  overcasts.  ©  2008  Optical  Society  of  America 
OCIS  codes:  010.1290,  010.1615,  010.1690,  010.7295,  290.1090,  330.1730. 


1.  Introduction 

Although  the  words  “overcast”  and  “color”  seem  un¬ 
likely  companions,  daytime  cloudy  skies  exhibit  col¬ 
ors  ranging  from  the  apparent  monochromes  of 
stratus  layers  to  the  remarkable  variegation  of  stra- 
tocumulus  (Sc)  overcasts.  Minnaert  conducted  some 
of  the  earliest  research  on  cloud  colors,  and  despite 
the  little  known  quantitatively  about  them  in  his 
day  [1],  he  nonetheless  offered  several  useful  quali¬ 
tative  insights.  Among  these  was  that  the  colors  of 
distant  cumuli  are  affected  by  the  changing  balance 
between  (1)  bluish  airlight  scattered  toward  an  ob¬ 
server  and  (2)  yellowish  extinction  of  scattered  cloud 
light  over  long  optical  paths.  Middleton  extended 
these  ideas  in  his  colorimetric  modeling  of  the  colors 
of  distant  objects  seen  through  clear  or  hazy  atmo¬ 
spheres  [2] ,  although  he  did  not  analyze  similar  col¬ 
ors  seen  beneath  overcasts.  This  early  work  has 
hardly  been  followed  by  a  flurry  of  activity:  papers 
that  measure  and  model  overcast  chromaticities  or 
visible-wavelength  spectra  are  rare  indeed  [3—6] , 
and  those  that  examine  the  angular  distribution  of 
overcast  colors  seem  nonexistent.  However,  recent 
advances  in  digital  image  analysis  now  let  us 
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measure  the  angular  color  variations  that  Minnaert 
could  only  outline. 

2.  Spectral  and  Colorimetric  Calibration  of  Skylight 
Cameras 

For  several  decades,  researchers  in  numerous  fields 
have  devised  techniques  for  recovering  quantitatively 
accurate  color  and  spectral  information  from  digitized 
photographs  [7—9] .  The  central  problem  faced  by  all 
such  techniques  is  calibrating  the  digital  imaging  sys¬ 
tem’s  response  to  a  wide  range  of  scene  spectral 
radiances  Lx  at  visible  wavelengths  L  In  trichromatic 
systems  such  as  red-green-blue  (RGB)  digital  still 
cameras,  the  recovered  spectra  are  only  metamers 
for  the  original  scene’s  LA.  This  is  true  because 
each  sensor  pixel’s  RGB  data  provide  only  three 
broadband  measures  of  the  spectral  power  distribu¬ 
tion,  rather  than  the  much  larger  number  of  LA  that 
typically  specify  the  original  (say,  30-100  LA  for  a  400- 
700  nm  visible  spectrum).  Strictly  speaking,  then,  all 
such  recovered  spectra  are  underdetermined. 

Yet  Nieves  et  al.  have  shown  recently  that  if  a  cali¬ 
brated  camera  photographs  only  scenes  whose  Lx  are 
similar  to  those  in  the  calibration  training  set,  then 
the  recovered  spectra  are  quite  realistic  both  colori¬ 
metrically  and  spectrally  [10].  Their  simple,  yet  com¬ 
putationally  robust,  calibration  algorithm  permits 
accurate  reconstructions  of  spectra  for  various 
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families  of  light  sources  (e.g.,  daylight  at  high  unre¬ 
fracted  sun  elevations  h0  or  skylight  at  low  h0).  The 
spectral  family  studied  here  is  overcasts  with 
h0  >  10°,  since  radiometer  spectra  acquired  during 
our  earlier  work  [6]  indicate  that  overcast  color 
changes  appreciably  at  lower  h0.  My  chief  modifica¬ 
tion  to  the  Nieves  calibration  procedure  is  to  build 
the  LA  training  set  from  directly  measured  rather 
than  diffusely  reflected  skylight  spectra. 

To  do  so,  I  use  a  Photo  Research  PR-650  spectrora- 
diometer  that  has  a  spectral  range  of  380-780  nm,  a 
step  size  of  4  nm  [11] ,  and  a  telescopic  lens  that  per¬ 
mits  radiance  measurements  across  a  1°  diameter 
field  of  view  (FOV).  As  shown  in  Fig.  1,  a  metal  bracket 
is  mounted  on  a  tripod  so  that  the  aims  of  the  attached 
cameras  and  radiometer  remain  fixed  with  respect 
to  one  another.  After  sighting  through  the  radio¬ 
meter’s  telescope  on  a  distant  target  and  then  photo¬ 
graphing  it,  I  can  match  the  radiometer’s  circular 
FOV  with  the  corresponding  1°  diameter  circle  on 
each  camera’s  sensor  array.  Then  the  radiometer’s 
LA  are  those  that  generate  the  camera’s  RGB  values 
averaged  across  the  1°  circular  FOV.  All  photographs 
are  taken  with  each  camera  recording  in  its 
12  bit/channel  RAW  mode  (which  has  a  fixed  spectral 
response),  and  the  only  image  processing  done  is  to 
multiply  all  RGB  values  by  a  constant  linear  gain  fac¬ 
tor.  This  fixed  stretching  of  image  histograms  makes 
photographs  legibly  bright  but  in  no  way  changes  the 
camera  spectral  response. 

Because  many  daylight  and  skylight  spectra  are 
fairly  simple  variants  on  extraterrestrial  sunlight, 
studies  consistently  show  that  the  first  few  eigenvec¬ 
tors  (typically  3-5)  of  a  principal-components  analy¬ 
sis  can  reproduce  these  spectra  with  great  accuracy 
[5,12,13].  The  calibration  algorithm  that  Nieves  et  al. 
call  their  “direct  pseudo-inverse  transform”  [10] 


Fig.  1.  Spectral  and  colorimetric  calibration  equipment  consists 
of  two  digital  still  cameras  (a  Nikon  E5000  and  an  Olympus  E300, 
labeled  “b”  and  “c,”  respectively)  and  a  Photo  Research  PR-650 
spectroradiometer  (labeled  “a”)  affixed  to  a  metal  bracket,  which 
is  itself  mounted  on  a  tripod. 


exploits  this  reduced  dimensionality  of  daylight 
and  skylight  spectra,  with  the  added  advantages 
that  one  need  not  (1)  calculate  eigenvectors  for  the 
training  set  of  spectra  or  (2)  determine  spectral 
sensitivities  for  a  camera’s  color  channels. 

Separate  tests  with  a  monochromator  and  a  nar¬ 
row  bandpass  filter  show  that  the  digital  cameras 
I  use  (see  Fig.  1)  respond  to  LA  throughout  the  visible, 
here  taken  to  be  400-700  nm.  At  the  radiometer’s 
4nm  resolution,  this  gives  76  per  visible  spectrum. 
To  reconstruct  such  spectra,  I  start  by  using  a  cam¬ 
era’s  three  unfiltered  RGB  channels  to  form  the  76  x 
3  transform  matrix  F  defined  by 

F  =  L  pT(ppT)-\  (1) 

where  L  is  a  76  x  n  matrix  of  n  different  skylight  ra¬ 
diance  spectra  measured  with  the  radiometer  (the 
training  set  of  spectra),  p  is  the  corresponding  3  xn 
matrix  of  normalized  12  bit/channel  RGB  pixel 
values  for  the  same  n  features,  and  pT  is  p’s  transpose. 
My  current  training  set  uses  measurements  from  five 
different  overcasts,  and  it  has  n  =  232  features  of 
widely  varying  color  and  brightness.  To  use  the  result¬ 
ing  least-squares  transform  matrix  F,  I  photograph  a 
new  overcast  feature  with  the  camera  still  operating 
in  its  native  RAW  mode.  This  new  feature’s  recon¬ 
structed  76  x  1  spectrum  Lr  is  then  given  by 

Lr  =  Fpr,  (2) 

where  pr  is  the  feature’s  3x1  matrix  of  RGB  pixel  va¬ 
lues.  Because  F  is  derived  from  similar  spectra,  the 
resulting  spectral  errors  in  Lr  will  be  small. 
Errors  in  CIE  1976  u\  v'  chromaticities  calculated 
from  Lr  can  be  reduced  ~12%-21%  by  applying  non¬ 
linear  spectral  corrections  to  Lr.  These  spectral  multi¬ 
pliers  are  calculated  by  (1)  forming  n  different  spectra 
of  the  wavelength-by-wavelength  ratios  of  L  :  Lr  in 
the  training  set  and  (2)  weighting  each  of  these  n  cor¬ 
rection  spectra  by  its  inverse  colorimetric  distance 
from  the  new  feature’s  Lr. 

The  high  quality  of  the  reconstructed  Lr  is  evident 
in  Fig.  2,  which  shows  the  95th  percentile  root-mean- 
square  error  (RMSE)  between  the  radiometer  and 
Nikon  E5000  camera  spectra  for  a  set  of  measure¬ 
ments  not  included  in  the  calibration  training  set. 
In  other  words,  only  5%  of  all  Lr  reconstructed  from 
Nikon  camera  data  can  be  expected  to  generate  a 
RMSE  larger  than  Fig.  2’s  small  0.22%  difference. 
Reconstructed  chromaticities  are  equally  accurate: 
when  compared  with  radiometer  spectra  not  in  the 
training  set,  the  corresponding  Nikon  Lr  yield  a 
mean  colorimetric  error  of  A u'v'  =  0.000359  and  a 
A u'v'  standard  deviation  of  0.000230  [14].  This  error 
is  clearly  subthreshold  because  it  is  <30%  of  the 
mean  MacAdam  just-noticeable  difference  (or  JND) 
[15]  for  all  overcast  chromaticities  measured  with 
the  radiometer. 
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Fig.  2.  Normalized  radiances  for  the  95th  percentile  RMSE  spec¬ 
trum  reconstructed  from  Nikon  E5000  camera  data  compared  with 
the  original  radiometer  spectrum  of  an  overcast  feature.  Each 
curve  is  normalized  so  that  its  spectral  radiances  Lx  satisfy 

E(^)2  =  1- 

3.  Colorimetry  of  Overcasts  from  All-Sky  Photographs 

Overcast  skies  yield  an  embarrassment  of  remote¬ 
sensing  riches:  cloud  optical  depths  t,  radiances, 
and  colors  all  can  vary  greatly  over  distances  of  a 
few  degrees  and  time  intervals  of  a  few  seconds. 
One  especially  useful  tool  for  capturing  meaningful 
details  from  this  endlessly  mutable  scene  is  a  combi¬ 
nation  of  fisheye  lens  and  camera  that  records  a  180° 
FOV.  After  equipping  Fig.  l’s  Nikon  E5000  with  a 
FC-E8  fisheye  converter  lens  [16],  my  midshipman 
research  students  and  I  at  the  U.  S.  Naval  Academy 
(USNA)  in  Annapolis,  Maryland,  took  more  than 
1300  such  fisheye  photographs  of  20  different  over¬ 
casts  from  August  2006  through  February  2008. 

As  in  our  earlier  work  [6,17] ,  overcasts  are  defined 
here  as  having  cloud  cover  sufficient  to  obscure  any 
visible  clear  sky  and  to  make  any  cast  shadows  indis¬ 
tinct.  Although  most  of  the  USNA  overcast  images 
were  of  nonprecipitating  Sc,  stratus,  and  altostratus, 
with  careful  lens  cleaning  it  was  also  possible  to 
photograph  overcasts  during  drizzle,  light  rain, 
and  snow.  The  resulting  all-sky  data  on  overcast 
spectra  and  colors  offer  several  surprises. 

A.  Azimuthally  Averaged  Chromaticity  Curves 

For  example,  Fig.  3  shows  a  closeup  view  of  the  CIE 
1976  UCS  diagram,  including  such  colorimetric  land¬ 
marks  as  part  of  the  Planckian  locus  and  its  corre¬ 
sponding  color-temperature  limits.  The  short 
horizontal  line  labeled  “1  JND”  is  a  perceptual  ruler 
whose  length  is  that  of  a  representative  MacAdam 
JND  for  Fig.  3’s  two  Sc  overcasts  at  USNA.  Each 
curve  extends  from  the  zenith  (view-elevation  angle 
h  =  90°)  to  near  the  horizon  (h  =  0°),  and  each  plots 
reconstructed  chromaticities  averaged  over  sky 


Fig.  3.  Portion  of  the  CIE  1976  UCS  diagram,  showing  Sc  over¬ 
cast  chromaticities  u',  v'  measured  at  USNA  in  Annapolis,  Mary¬ 
land,  on  30  August  2006  and  12  January  2007.  Each  curve  plots 
temporally  and  azimuthally  averaged  chromaticities  versus 
view-elevation  angle  h  above  the  astronomical  horizon;  the  mean 
unrefracted  sun  elevation  h0  during  each  curve’s  photography 
session  (here,  ~38  min)  is  also  listed. 

azimuths  (i.e.,  along  almucantars).  Because  overcast 
radiances  are  often  taken  to  be  uniform  along 
almucantars  [18],  this  choice  of  averaging  direction 
is  reasonable  if  not  unassailable  [19]. 

To  reduce  radiance  and  color  fluctuations  at  each 
pixel  caused  by  short-term  changes  in  cloud  t  [17], 
photographs  were  taken  at  30  s  intervals  and  aver¬ 
aged  over  ^38  min  in  Fig.  3  to  form  a  mean  image 
of  each  overcast  (see  Fig.  4).  In  my  experience,  such 
time  periods  are  long  enough  to  produce  useful  colori¬ 
metric  smoothing  yet  brief  enough  to  capture  many 
important  short-term  overcast  features.  Two  poten¬ 
tially  misleading  details  are  noteworthy  in  Fig.  4. 
First,  averaging  makes  these  mean  Sc  images  resem¬ 
ble  stratus.  Second,  the  faint  arcs  that  stretch  from 
horizon  to  horizon  in  Fig.  4(b)  are  not  crepuscular 
rays  but  instead  are  caused  by  persistent  maxima 
in  cloud  t  along  the  mean  flow,  and  these  arcs  are 
sometimes  accompanied  by  gravity  waves  nearly  per¬ 
pendicular  to  them. 

Each  Fig.  3  chromaticity  curve  derived  from  Fig.  4 
has  a  pair  of  extrema,  with  one  such  colorimetric 
turning  occurring  at/i~33°or/i~380  and  the  other 
at  h  ~  2.25°  or  h  ~  2.75°,  respectively.  While  these  h 
values  vary  widely  across  overcasts,  the  turnings 
themselves  are  quite  common  features:  as  we  look 
from  high  to  low  h,  the  first  turning  is  counterclock¬ 
wise,  and  the  second  is  clockwise.  For  the  12  January 
2007  overcast,  Fig.  3  restricts  azimuthal  averages  to 
a  167°  wide  sector  that  includes  a  bright  yellowish 
overcast  area  near  the  horizon  (see  the  10  o’clock  po¬ 
sition  in  Fig.  4(b)),  and  this  area  accounts  for  Fig.  3’s 
comparatively  large  color  gamut  at  low  h.  Obviously 
this  brightening  is  caused  by  reduced  cloud  t  or  even 
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Fig.  4.  Time-averaged  fisheye  images  of  nonprecipitating  Sc  overcasts  photographed  at  USNA  on  (a)  30  August  2006  and  (b)  12  January 
2007.  The  mean  image  in  (a)  is  calculated  from  53  individual  photographs  and  has  a  measured  cloud  base  z  ~  365  m,  whereas  (b)  is  cal¬ 
culated  from  78  photographs  and  has  2400  m  <  z  <  2700  m. 


by  breaks  in  the  overcast  in  these  directions,  but  so 
long  as  no  clear  sky  is  visible  anywhere  the  cloud  cov¬ 
er  still  satisfies  my  definition  of  overcast.  At  least  at 
USNA,  such  horizon  brightening  is  not  rare  in  Sc 
overcasts,  and  it  occurs  over  both  land  and  the 
adjacent  Chesapeake  Bay 
Like  Fig.  3,  Fig.  5  also  plots  meridional  curves  of 
overcast  chromaticities  that  are  averaged  temporally 
and  azimuthally.  Note  that  although  both  Figs.  3  and 
5  are  scaled  isotropically  in  u'  and  v' ,  Fig.  5  spans  a 
much  smaller  colorimetric  area  (A u!  =  Auf  =  0.015 
in  Fig.  5  as  opposed  to  0.025  in  Fig.  3).  Thus,  while 
curvatures  can  be  compared  directly  in  these  figures, 
Table  1  shows  that  their  normalized  chromaticity 
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Fig.  5.  Temporally  and  azimuthally  averaged  chromaticity 
curves  for  Sc  overcasts  measured  at  USNA  on  20  February  and 
7  March  2007. 


gamuts  g  are  distinctly  different  [20].  Table  l’s 
conventional  meteorological  abbreviations  “OVC” 
and  “CLR”  indicate  overcast  and  clear,  respectively. 

Differences  in  g  aside,  nearly  every  Sc  overcast 
that  we  photographed  has  two  distinct  colorimetric 
turnings,  one  at  high  h  and  the  other  at  low  h.  How¬ 
ever,  the  turnings’  details  can  vary  considerably 
from  day  to  day,  just  as  they  do  in  the  clear  sky. 
For  example,  Fig.  6  pairs  an  overcast  with  large, 
nearly  symmetric  turnings  (9  November  2007  curve) 
with  one  where  these  features  are  vestigial  (29  No¬ 
vember  2006  curve).  A  few  overcasts  are  bluest  near 
the  horizon,  and  this  causes  particularly  large  colori¬ 
metric  excursions  toward  the  blue  at  low  h.  The  30 
August  2006  and  7  March  2007  Sc  overcasts  (Figs.  3 
and  5,  respectively)  clearly  show  this  fairly  unusual 
behavior. 

A  simple  semiempirical  model  suggests  that  this 
low-/i  turning  toward  the  blue  (i.e.,  toward  higher  col¬ 
or  temperatures;  see  Section  4)  is  caused  by  bluish 
airlight  whose  radiances  steadily  increase  with  de¬ 
creasing  h  until  just  above  the  horizon  [21].  Perhaps 
surprisingly,  overcasts  thus  repeat  the  clear  sky’s 
pattern  of  skylight  reddening  at  intermediate  h ,  fol¬ 
lowed  by  airlight  bluing  at  lower  h.  Figure  7  illus¬ 
trates  this  colorimetric  behavior  common  to  clear 
and  overcast  skies,  although  here  the  turning’s  direc¬ 
tion  at  high  h  does  differ  for  the  two  sky  states  (i.e., 
clockwise  with  decreasing  h  in  the  clear  sky  and 
counterclockwise  in  the  overcast).  Such  subtle  trends 
would  be  difficult  indeed  to  observe  without  the  de¬ 
tailed  data  provided  by  all-sky  digital  imaging. 

B.  Histograms  of  Chromaticity  and  Inverse  Color 
Temperature 

We  can  see  even  more  colorimetric  detail  in  overcasts 
by  plotting  two-dimensional  histograms  of  their 
chromaticities.  Although  this  added  detail  can  be 
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Table  1.  Mean  Skylight  Chromaticities  and  Chromaticity  Gamuts 


Figure 

Date 

Sky  State 

h0  Interval6  (°) 

Mean  u' 

Mean  v' 

8 

3 

8-30-06 

OVC 

46.4-52.1 

0.19883 

0.47008 

0.0034196 

3“ 

1-12-07 

OVC 

27.6-29.1 

0.20037 

0.47379 

0.02275 

5 

2-20-07 

OVC 

40.1-39.9 

0.1995 

0.47063 

0.012948 

5a 

3-7-07 

OVC 

39.6-43.2 

0.19865 

0.46834 

0.0042835 

6 

11-29-06 

OVC 

25.2-27.5 

0.19954 

0.47145 

0.0053575 

6 

11-9-07 

OVC 

15.6-10.1 

0.1964 

0.4636 

0.0068814 

7 

1-17-07 

CLR 

29.8 

0.1783 

0.42476 

0.065473 

7 

10-5-06 

OVC 

15.2-9.4 

0.19641 

0.4665 

0.0073714 

8  a 

1-12-07 

OVC 

27.6-29.1 

0.19883 

0.4706 

0.014246 

9a 

3-7-07 

OVC 

39.6-43.2 

0.1986 

0.46806 

0.0058461 

11 

9-14-07 

OVC 

53.0-50.6 

0.19905 

0.47053 

0.0069569 

aFor  a  given  overcast,  gamut  g  and  mean  u' ,  v'  depend  slightly  on  whether  chromaticities  are  first  averaged  azimuthally  (Figs.  3  and  5)  or 
not  (Figs.  8  and  9). 

^Morning  observations  list  the  smaller  h0  value  first,  whereas  afternoon  observations  list  it  last. 


instructive,  it  comes  at  the  price  of  obscuring  the 
angular  distributions  of  overcast  colors.  For  example, 
Fig.  8  is  a  histogram  of  all  colors  in  the  12  January 
2007  overcast  (i.e.,  sans  azimuthal  averaging;  com¬ 
pare  Fig.  3),  and  its  gray-scale  legend  scales  relative 
frequency  fre[  logarithmically.  Thus  the  bright  pixels 
clustered  near  u'  =  0.198,  v'  =  0.469  indicate  the  pre¬ 
dominance  of  those  colors  in  this  overcast.  However, 
equally  evident  is  Fig.  8’s  curve  of  low-/^  pixels  that 
arcs  away  from  the  Planckian  locus  and  ends  in 
Fig.  4(b)’s  small  but  distinct  band  of  yellowish 
overcast  near  the  horizon. 

Figure  8  provides  information  not  found  in  Fig.  3, 
such  as  the  fact  that  overcast  chromaticities  can  be 
on  the  purplish  or  right-hand  side  of  the  Planckian 
locus.  Not  surprisingly,  Fig.  8’s  individual  chromati¬ 
cities  span  a  larger  area  on  the  CIE  diagram  than  do 
their  azimuthally  averaged  kin  in  Fig.  3.  Neither  fig- 


Fig.  6.  Temporally  and  azimuthally  averaged  chromaticity 
curves  for  Sc  overcasts  measured  at  USNA  on  29  November 
2006  and  9  November  2007.  Light  drizzle  fell  throughout  measure¬ 
ments  on  the  latter  date.  To  show  more  detail,  the  u' ,  v'  scaling  is 
slightly  anisotropic  here. 


ure  likely  jibes  with  our  visual  impression  of  Fig.  4(b) 
— isolated  bands  of  yellow  near  the  horizon  appear  to 
be  paired  with  near  monochromes  elsewhere.  Yet 
Figs.  3  and  8  together  make  clear  that  even  this 
seeming  dichotomy  of  color  is  in  fact  a  continuum. 

As  Table  1  indicates,  Fig.  8’s  overcast  chromaticity 
gamut  g  is  unusually  large.  More  typical  is  Fig.  9’s 
modest  array  of  chromaticities  for  the  much  more 
uniform  Sc  overcast  of  7  March  2007,  for  which  g 
is  only  ~40%  of  its  value  on  12  January  2007.  This 
reduced^  stems  from  the  absence  in  Fig.  9  of  Fig.  8’s 
near-horizon  yellowing,  which  is  itself  caused  by 
small  areas  of  thinner  overcast  in  Fig.  4(b).  Even 
though^  is  much  smaller  in  Fig.  9,  the  accompanying 
JND  scale  indicates  that  differences  are  still  discern¬ 
ible  among  its  overcast  colors. 


CIE  1976  u' 

Fig.  7.  Temporally  and  azimuthally  averaged  chromaticity 
curves  for  a  clear  (CLR)  sky  and  a  Sc  overcast  (OVC)  measured 
at  USNA  on  17  January  2007  and  5  October  2006,  respectively. 
The  CLR  chromaticities  are  averaged  within  two  90°  wide  sectors 
that  are  symmetric  about  the  clear-sky  principal  plane  (relative 
azimuths  of  45°-135°  and  225°-315°),  so  these  colors  exclude 
the  solar  and  antisolar  regions.  To  show  as  much  detail  as  possible, 
the  u',  v'  scaling  is  anisotropic  here. 
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Fig.  8.  Histogram  of  u',  v'  chromaticities  for  Sc  overcast  mea¬ 
sured  at  USNA  on  the  morning  of  12  January  2007  when 
27.6°  <  h0  <  29.1°.  The  absolute  pixel  frequency  in  any  u',  v1  bin 
is  the  product  of  its  gray-scale  relative  frequency  frel  and  the  total 
number  of  image  pixels  n. 


Note  that  the  modal  colors  in  Figs.  8  and  9  occur 
near  the  blue  ends  of  their  chromaticity  clusters, 
meaning  that  outlier  (i.e.,  low-frequency)  chromatici¬ 
ties  in  their  overcasts  tend  to  be  yellower  rather  than 
bluer.  However,  Figs.  10  and  11  show  that  the  reverse 
can  be  true:  in  the  14  September  2007  overcast,  out¬ 
lier  chromaticities  are  bluer  rather  than  yellower,  as 
indicated  by  Fig.  ll’s  swath  of  bright  pixels  to  the  left 
of  its  modal  chromaticities  at  u'  =  0.199,  v'  =  0.471. 
These  outlier  chromaticities  occur  mostly  in  Fig.  IP's 
first  trigonometric  quadrant  just  above  its  horizon. 
Other  bluish  areas  in  Fig.  10  include  the  thin  arcs 


Fig.  9.  Histogram  of  u',  v'  chromaticities  for  Sc  overcast  mea¬ 
sured  at  USNA  on  the  morning  of  7  March  2007  when 
39.6°  <  h0  <  43.2°.  Throughout  measurements,  precipitation  var¬ 
ied  from  light  snow  to  flurries. 


Fig.  10.  Time-averaged  fisheye  image  calculated  from  61  indivi¬ 
dual  photographs  of  a  Sc  overcast  photographed  at  USNA  on  14 
September  2007;  the  measured  cloud  base  is  ~2400  m.  For  consis¬ 
tency  in  printing,  no  image  processing  is  used  to  correct  the  small, 
but  constant,  color  bias  of  this  and  other  RAW-format  images 
(Figs.  4(a)  and  4(b)). 


that  stretch  across  the  sky  and  are  slightly  darker 
than  their  surroundings.  As  is  true  in  Fig.  4(b),  the 
simplest  explanation  for  these  long-lived  arcs  is  per¬ 
sistent  maxima  in  cloud  t  along  the  mean  flow.  This 
increased  cloud  t  both  decreases  the  arcs'  radiances 
and  increases  their  color  temperatures  (i.e.,  makes 
them  slightly  bluer)  [6] . 

Correlated  color  temperature  (CCT)  has  long 
served  as  an  approximate,  convenient  alternative 
to  chromaticity  coordinates  near  the  Planckian  locus 
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Fig.  11.  Histogram  of  u',  v'  chromaticities  for  nonprecipitating  Sc 
overcast  measured  at  USNA  on  the  afternoon  of  14  September 
2007  when  53.0°  >  h0  >  50.6°. 
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Hill 


[22].  Because  CCT  reduces  two-dimensional  u! ,  1/ 
data  to  one  dimension,  CCT  necessarily  introduces 
some  colorimetric  ambiguity  into  quantifying  over¬ 
cast  color  [23].  For  example,  an  overcast  feature 
whose  CCT  =  6000  K  can  lie  on  either  side  of  the 
Planckian  locus.  But  since  the  perceptual  conse¬ 
quences  of  this  uncertainty  often  are  moot  for  over¬ 
cast  colors,  CCTs  are  useful  here  provided  that  we 
keep  their  limitations  in  mind.  Even  more  useful 
than  CCT  itself  is  inverse  CCT,  which  yields  a  scale 
that  describes  human  color  response  more  uniformly 
[5].  The  unit  of  inverse  CCT  is  the  reciprocal  mega- 
Kelvin  (106/CCT),  which  is  denoted  by  the  symbol 
MK"1  and  was  originally  called  the  “mired”  [24]. 

Figure  12  interleaves  histograms  of  inverse  CCT 
for  the  7  March  2007  and  14  September  2007  Sc  over¬ 
casts;  note  that  the  ordinate  is  log10(n),  where  n  is 
the  number  of  pixels  in  each  inverse-CCT  interval. 
On  comparing  Fig.  12’s  inverse-CCT  histograms  with 
the  corresponding  chromaticity  histograms  in  Figs.  9 
and  11,  notice  how  the  bluer  colors  (i.e.,  those  with 
smaller  inverse  CCTs)  occur  with  greater  relative 
frequency  on  14  September  2007  than  they  do  on  7 
March  2007.  In  other  words,  the  14  September 
2007  histogram  in  Fig.  12  is  more  symmetric  about 
its  modal  color  than  is  the  7  March  2007  histogram. 
Both  the  bluish  horizon  and  faint  arcs  in  Fig.  10  con¬ 
tribute  to  this  enhanced  blueness. 

In  Fig.  13,  the  abscissa  scale  changes  in  order  to 
interleave  the  7  March  2007  and  12  January  2007 
histograms,  a  comparison  that  clearly  reveals  the  lat¬ 
ter  overcast’s  marked  color  skewness.  Figure  13’s 
long  tail  of  overcast  colors  with  lower  CCTs  (and  thus 
higher  inverse  CCTs)  on  12  January  2007  is  just  as 
prominent  and  easy  to  quantify  here  as  it  is  in  Fig.  8. 
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Fig.  12.  Histograms  of  inverse  CCT  for  overcasts  measured  at 
USNA  on  7  March  2007  (solid  bars)  and  14  September  2007 
(striped  bars).  Each  bar’s  height  is  proportional  to  log10(ft),  where 
n  is  the  number  of  pixels  in  the  corresponding  interval  of  inverse 
CCT.  Each  juxtaposed  pair  of  bars  for  the  two  dates  has  a  common 
abscissa  value  (e.g.,  the  leftmost  pair  is  for  144  MK-1). 


Fig.  13.  Histograms  of  inverse  CCT  for  overcasts  measured  at 
USNA  on  7  March  2007  (solid  bars)  and  12  January  2007  (striped 
bars). 

Once  again,  though,  graphical  clarity  should  not  be 
confused  with  easy  visibility:  the  continuum  of  colors 
present  in  Fig.  4(b)  remains  as  difficult  to  quickly  dis¬ 
cern  there  as  it  is  immediately  obvious  in  Figs.  8  and 
13.  Yet  one’s  real-time  experience  of  a  constantly 
changing  Sc  overcast  need  not  be  divorced  from  such 
colorimetric  order.  Indeed,  my  experience  is  that 
foreknowledge  makes  such  visual  subtleties  easier 
to  see. 

4.  Modeling  the  Angular  Distribution  of  Overcast 
Colors 

Our  past  experience  with  general-purpose  radiative 
transfer  models  such  as  MODTRAN4  suggests  that 
they  can  simulate  individual  spectra  of  overcast  ra¬ 
diances  fairly  realistically.  This  qualified  endorse¬ 
ment  means  that,  for  a  reasonable  choice  of  input 
parameters,  MODTRAN4  produces  spectra  that 
are  nearly  congruent  with  individual  measured  over¬ 
cast  spectra  [6].  However,  neither  MODTRAN4  nor 
its  predecessor  LOWTRAN7  appears  to  simulate 
the  angular  patterns  of  overcast  spectra  and  chroma- 
ticities  particularly  well. 

Consider  the  divergence  between  these  models  and 
reality  in  Fig.  14.  For  the  sake  of  comparison,  Fig.  14 
repeats  Fig.  5’s  meridional  curve  of  the  azimuthally 
averaged  colors  that  were  measured  in  the  20  Febru¬ 
ary  2007  Sc  overcast.  Paired  with  these  observations 
are  comparable  overcast  chromaticity  curves  from 
LOWTRAN  and  MODTRAN.  Each  simulation  uses 
its  default  Sc  cloud  model,  with  cloud  bases  set  to 
either  z  =  0.66  km  or  z  =  1.0  km  above  the  surface 
and  cloud  thicknesses  A z  =  1.34  km  or  A z  =  1.9  km 
(LOWTRAN  and  MODTRAN  parameters,  respec¬ 
tively).  For  each  model,  h0  =  45°,  and  the  viewing 
azimuth  is  90°  from  the  sun.  Using  only  a  single  re¬ 
lative  azimuth  is  acceptable  because  changing  this 
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curve  for  the  20  February  2007  Sc  overcast  compared  with  Sc  si¬ 
mulations  by  LOWTRAN7  and  MODTRAN4.  For  a  wide  range  of 
input  parameters,  neither  model  predicts  overcast  chromaticity 
patterns  well  and  both  greatly  underestimate  typical  overcast 
color  gamuts. 


parameter  affects  the  models’  chromaticities  negligi¬ 
bly  for  such  optically  thick  overcasts. 

These  models’  overcast  chromaticity  patterns  are 
odd  indeed.  Although  all  the  overcasts  that  we 
measured  exhibit  just  the  opposite  color  trend, 
LOWTRAN  predicts  consistent  reddening  of  the  over¬ 
cast  from  the  zenith  down  to  h~  6°,  and  only  there 
does  bluing  begin.  No  more  satisfactory  is 
MODTRAN’s  claim  that  the  Sc  overcast  becomes  pro¬ 
gressively  bluer  from  the  zenith  until  h  ~  0.5°,  with  a 
minuscule  move  toward  the  red  thereafter.  At  such 
low  h ,  MODTRAN’s  assumption  of  a  plane-parallel  at¬ 
mosphere  is  problematic  at  best  and  may  well  contri¬ 
bute  to  its  nonphysical  behavior.  Note  too  that 
MODTRAN’s  overcast  color  gamut  is  much  smaller 
than  either  those  predicted  by  LOWTRAN  or  ob¬ 
served  outdoors.  Changing  overcast  z  and  A z  does 
not  affect  the  shape  of  either  model’s  chromaticity 
curves  much,  which  suggests  that  their  verisimilitude 
problems  cannot  be  solved  simply  by  adjusting  these 
parameters.  Given  the  complexity  of  these  models, 
the  causes  of  their  colorimetric  problems  are  not 
immediately  obvious. 

Recall  that  observed  overcast  color  patterns  typi¬ 
cally  consist  of  (a)  bluing  from  the  zenith  down  to 
h  ~  30°-40°,  followed  by  (b)  reddening  down  to 
h  ~  3°,  and  finally  (c)  bluing  down  to  the  horizon. 
Figure  15  illustrates  these  three  distinct  color  re¬ 
gimes  as  seen  in  overcasts  photographed  at  mark¬ 
edly  different  h0  on  5  October  2006  and  20 
February  2007.  A  simple,  physically  plausible  expla¬ 
nation  for  these  regimes  is  that  surface-based  obser¬ 
vers  see  a  shifting  balance  between  the  spectral 
effects  of  (1)  absorption  within  the  overcast  and  both 


and  the  changing  balance  of  radiative  transfer  factors  that  ex¬ 
plains  them:  (1)  bluing  caused  by  increasing  spectral  absorption 
in  clouds  from  the  zenith  to  h  ~  30°-40°,  (2)  reddening  caused 
by  increasing  extinction  beneath  the  overcast  over  ever-longer 
slant  optical  paths  for  30°-40°  >  h  >  3°,  and  (3)  bluing  caused 
by  increasing  airlight  for  h  <  3°. 


(2)  extinction  losses  and  (3)  airlight  gains  that  arise 
from  scattering  in  the  cloud-free  air  beneath  it.  Pro¬ 
cess  (1)  causes  the  initial  bluing  from  the  zenith  to 
h  ~  30°-40°,  while  the  latter  two  yield  (2)  reddening 
until  a  few  degrees  above  horizon  and  then  (3)  bluing 
at  still  lower  h. 

My  semiempirical  model  attempts  to  describe  over¬ 
cast  colors  realistically  by  combining  spectrally  selec¬ 
tive  absorption  and  scattering  in  the  overcast  and 
subcloud  layers.  Its  theoretical  basis  is  the  radiative 
transfer  equation  [25],  here  given  as 

Lf  =  L0  exp(-T^)  +  /  f  (J(t)  exp(r  -  Tf))dr,  (3) 
Jo 

A  B 

where  the  radiance  Lf  received  by  an  observer  de¬ 
pends  on  the  combined  effects  of  terms  A  and  B.  Term 
A  is  the  attenuated  radiance  L0  for  a  distant  light 
source  (here  the  cloud  base;  see  Fig.  16),  and  term 
B  is  an  airlight  term  that  represents  the  cumulative 
scattering  effects  of  differential  atmospheric  volumes 
dV  along  any  slant  path  from  t  =  0  at  L0  to  ty,  the 
total  optical  depth  along  the  path  between  observer 
and  cloud  base.  A  source  function  J (t)  describes  the 
scattering- angle  dependence  of  dV’s  illumination  by 
skylight  radiances  and  the  net  radiance  that  dV 
scatters  toward  the  observer.  All  terms  in  Eq.  (3) 
are  implicit  functions  of  wavelength  A. 

At  visible  A  below  the  cloud  base,  term  A  tallies 
source  radiances  not  lost  due  to  extinction  by  scatter¬ 
ing  along  Tf,  while  term  B  tallies  gains  of  scattered  ra¬ 
diance  from  the  intervening  air.  Few  real  atmospheric 
J(r)  yield  airlight  integrals  (term  B)  that  can  be 
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Fig.  16.  Scattering  geometry  for  an  object  of  radiance  L0  seen  by 
an  observer  through  a  multiple-scattering  medium  along  a  path  of 
total  optical  depth  t f.  The  path  consists  of  differential  air  volumes 
dV,  each  having  an  angular  scattering  phase  function  J(t).  Here 
L0  is  the  spectral  radiance  emerging  from  the  cloud  base  and  Tf  is 
for  the  nominally  clear  air  beneath  it. 

evaluated  analytically.  But  in  the  shallow  air  layer  be¬ 
neath  many  overcasts,  we  can  reasonably  assume 
that  the  J (t)  spectrum  does  not  change  significantly 
along  a  given  slant  path.  In  that  case,  the  spectral 
radiative  transfer  equation  reduces  to  LfX  =  L0A 
exp(-Ty^)  +  JA{  1  -  exp(-r^)). 

The  overcast  model’s  empirical  bases  are  mea¬ 
sured  values  of  (1)  the  cloud  spectral  transmissivities 
T ;l  that  yield  L0  [6],  (2)  the  subcloud  spectral  extinc¬ 
tion  coefficients  (3exth  with  Tf  =  /?exts  and  s  =  the 
slant-path  geometrical  distance  to  the  cloud  base, 
and  (3)  the  subcloud  J \  spectrum,  which  varies  with 
h  but  not  along  s  at  a  given  h. 

What  insights  does  this  model  yield?  Figure  17 
shows  that,  like  the  measured  spectra,  it  produces 
three  distinct  chromaticity  regimes.  For  the  purposes 
of  this  initial  demonstration,  the  model’s  chromati¬ 
city  trends  are  more  important  than  its  absolute 
chromaticities:  subsequent  refinements  in  its  empiri¬ 
cal  parameters  will  make  the  observed  and  simu¬ 
lated  curves  more  nearly  congruent.  In  regime  (a), 
from  the  zenith  to  h~  30°-40°,  the  model  overcast 


Fig.  17.  Temporally  and  azimuthally  averaged  chromaticity 
curve  for  the  20  February  2007  Sc  overcast  compared  with  a  simu¬ 
lation  from  the  semiempirical  model.  Like  MODTRAN  and  LOW- 
TRAN,  this  model  is  based  on  the  radiative  transfer  equation  and, 
unlike  them,  it  correctly  reproduces  the  chromaticity  gamuts  and 
turnings  observed  in  real  overcasts. 


becomes  progressively  bluer  with  decreasing  h  be¬ 
cause  here  the  cloud  absorption  spectrum  dominates 
and  increases  the  bluish  bias  of  L0.  This  L0  trend 
holds  at  all  elevation  angles  because  geometrical 
and  optical  slant  paths  through  the  overcast  steadily 
increase  as  h  decreases.  But  only  in  regime  (a)  does 
spectrally  selective  absorption  by  cloud  droplets  gov¬ 
ern  the  overcast’s  color  as  seen  from  the  surface. 

In  regime  (b),  for  30°-40°  >  h  >  3°,  net  reddening 
occurs  because  here  the  subcloud  extinction  spectrum 
[the  exp(-Tf)  factor  in  Eq.  (3)’s  term  A]  dominates. 
This  is  analogous  to  the  progressive  yellowing  with 
distance  of  cumulus  clouds  on  a  clear  day  [1,2,26], 
especially  when  the  sun  is  low  in  the  sky:  as  subcloud 
Tf  increases  with  decreasing  h ,  we  can  begin  to  see  the 
consequences  of  clear-air  scattering  preferentially  re¬ 
moving  shorter  wavelengths  from  the  unscattered  L0 
reaching  us.  Finally,  progressive  bluing  occurs  for  h  < 
3°  in  regime  (c):  at  these  largest  subcloud  Tf ,  ever- 
increasing  multiple  scattering  results  in  the  predomi¬ 
nance  of  bluish  airlight  [Eq.  (3)’s  term  B].  Taken 
together,  the  model’s  absorption  and  scattering  re¬ 
gimes  yield  a  perhaps  unexpected  result:  overcast 
and  clear  skies  have  as  many  radiative  transfer 
similarities  as  they  do  differences. 

5.  Conclusions 

These  research  results  address  two  long-standing 
problems  in  the  optics  of  overcast  skies:  capturing 
their  surfeit  of  angular  and  temporal  details  while 
trying  to  divine  some  order  from  their  seeming  chaos. 
Our  earlier  time-series  analyses  of  overcast  color  [6] 
and  brightness  [17]  stressed  overcasts’  temporal 
patterns,  while  here  the  emphasis  is  on  their  time- 
averaged  angular  distributions  of  color.  Naturally 
the  chromaticity  gamuts  and  details  of  these  colors 
change  perceptibly  from  one  overcast  to  the  next. 
Yet  their  underlying  meridional  structure  is  a  re¬ 
markably  consistent  set  of  three  distinct  chromati¬ 
city  regimes:  (a)  bluing  from  the  zenith  down  to 
h  ~  30°-40°,  then  (b)  reddening  down  to  h  ~  3°, 
and  finally  (c)  bluing  down  to  the  horizon.  Further¬ 
more,  the  colorimetric  turnings  that  separate  these 
overcast  regimes  are  nearly  always  the  same:  for  de¬ 
creasing  h ,  the  first  turning  is  counterclockwise  and 
the  second  is  clockwise.  Although  clear  skies  exhibit 
the  same  kind  of  turnings  (Fig.  7),  the  one  at  higher  h 
often  is  opposite  its  overcast  counterpart. 

Unlike  existing  models,  the  semiempirical  model 
presented  above  helps  explain  this  and  other  over¬ 
cast  color  features.  In  both  overcast  and  clear  skies, 
the  colorimetric  turnings  near  the  horizon  seem  to 
result  from  a  shift  between  the  predominance  of 
reddish  extinction  at  higher  h  and  bluish  airlight 
at  lower  h.  In  all  such  color  regimes,  scattering  by 
cloud-free  air  determines  what  we  see.  Yet  the 
high-h  turnings  differ  for  the  two  sky  states  because 
only  in  the  overcast  do  we  see  a  transition  from  the 
effects  of  cloud  absorption  to  those  of  clear-air  scat¬ 
tering.  Such  connections  between  the  optics  of  clear 
and  overcast  skies  appear  obvious  in  hindsight,  but 
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they  would  be  difficult  indeed  to  make  without  the 
new,  detailed  views  of  overcasts  made  possible  by 
all-sky  imaging. 
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